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The motion of a single hole in a Mott antiferromagnet is investigated based on the t ~ J model. 
An exact expression of the energy spectrum is obtained, in which the irreparable phase string effect 
[Phys. Rev. Lett. 77, 5102 (1996)] is explicitly present. By identifying the phase string effect with 
' spin backflow, we point out that spin-charge separation must exist in such a system: the doped 

hole has to decay into a neutral spinon and a spinless holon, together with the phase string. We 
' show that while the spinon remains coherent, the holon motion is deterred by the phase string, 

I resulting in its localization in space. We calculate the electron spectral function which explains the 

line shape of the spectral function as well as the "quasiparticle" spectrum observed in angle-resolved 
photoemission experiments. Other analytic and numerical approaches are discussed based on the 
present framework. 
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I. INTRODUCTION 

Since the discovery of high- Tc superconductivity, there has been much effort to elucidate the properties of doped 
Mott insulators. In this context, the specific case of one hole in a Mott insulating antiferromagnet (half filled band) has 



I 

(— I , been the subject of extensive theoretical and experimental investigation. These studies address the basic question, 
Q ■ whether the motion of a hole in the antiferromagnet can possibly be described within a quasiparticle approach. 
O , Photoemission spectroscopy provides valuable information that can help resolve this issue. Experimental results from 
angle resolved photoemission spectroscopy (ARPES) are now available for Sr2Cu02Cl2 as well as for Ca2Cu02Cl2 
: 10^- Both these materials are Mott insulators and are parent compounds of the high-Tc cuprates. The results from 
K*" ARPES can be summarized as follows: (i) the spectral features observed are not sharp at all. Quite to the contrary, an 
intrinsic broad feature extending to energies of the order of 1.5 eV before merging into the main valence band is seen. 
This is to be contrasted with sharp spectral features that one expects in a quasiparticle scenario; (ii) the observed 
, dispersion is isotropic around kg = (7r/2,7r/2); (iii) the measurements of Ronning et al. reveal the presence of a 
■ so-called remnant Fermi surface of the momentum structure in the energy-integrated spectral function. 
' The broad spectral features seen in ARPES strongly suggest the breakdown of the quasiparticle picture. Based 
on such considerations, Laughlin Q conjectured the failure of quasiparticle theory in these materials, and further 
, proposed, that the observed isotropic dispersion has its origins in an underlying spinon spectrum. This picture 

a envisages spin-charge separation with the scale of the observed dispersion determined by the superexchange J. 
I An entirely different picture is presented by the self-consistent Born approximation (SCBA) approach 

Though this scheme is based on the t — J model, it depicts a spin-polaron picture for the single hole case where 
^ \ the doped hole behaves not very different from a quasiparticle in the Landau-Fermi liquid theory. The SCBA results 
O ■ have also been supported by exact numerical calculations on finite lattices [0 up to 32 sites as well as variational 
calculations on larger lattices IT^ . Here, the consistency amongst different numerical methods mostly concerns the 
spectrum ek, usually defined as the minimum energy for a given momentum k. The spectrum is anisotropic around 
the "Fermi points" kg = (±7r/2, ±7r/2) in all these calculations. It is called a "quasiparticle" spectrum since a sharp 
; ^ ■ peak in the spectral function usually appears at ek in both the SCBA as well as results from exact diagonalization. 
5^ \ Note that these results are not consistent with ARPES which exhibits an isotropic dispersion around ko, and rather 
broad spectral features. Thus, to account for the former, the inclusion of second (t') and third {t") nearest neighbor 
hoppings to the t — J model have been proposed in the literature . However, this approach is meaningful 

only when the quasiparticle description holds, i.e., only if a sharp quasiparticle peak exists at ek in the t — J model. 
As illustrated in Fig. 1, if in the thermodynamic limit, the spectral function for a given momentum k does not show 
a sharp quasiparticle peak, then, notwithstanding how accurately ek is determined by various theoretical methods, 
it has no experimental implication. For, in this case, the higher energy e[. at the broad "peak" in Fig. 1, which is 
observable experimentally, may have nothing to do with the anisotropic ek. 

On more general grounds, the breakdown of the Landau- Fermi quasiparticle picture has been discussed by Anderson 
for a doped Mott insulator in the presence of an upper Hubbard band [0. He argues that the quasiparticle weight 
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FIG. 1. Schematic illustration of a spectral function at a given k in which the quasiparticle weight ~ a.t the energy 
bottom ek. Here the "quasiparticle" peak at ej^ h^-s nothing to do with the Landau- Fermi quasiparticle. 



Z vanishes due to unrenormalizable Fermi-surface phase shifts induced when a particle is injected into the Mott 
insulator. Indeed, based on a rigorous formulation of the single-electron Green's function in the one-hole case using 
the t — J model, it has been demonstrated that the quasiparticle weight Z at the Fermi surface must vanish in 

the thermodynamic limit due to the presence of the phase string effect. Such a phase string effect can be considered 
as the equivalent of the phase shifts proposed by Anderson in the Hubbard model. 

This paper concerns the effect of the aforementioned phase string on the dynamics of a single hole introduced in 
an antiferromagnet. We show how the phase string leads to the frustration of kinetic energy of the hole. We then 
show that the phase string effect is related to spin-charge separation and that the charge carrier is actually a spinless 
holon. We find that the holon propagator is localized in space, owing to the phase string. This is in contrast to the 
prediction of SCBA that the hole behaves as a Landau quasiparticle. Our result for the holon propagator is consistent 
with an earlier conclusion that the quasiparticle weight is zero for the doped hole ||l^,|l9|. It does not necessarily 
contradict finite-size calculations, as the localization length scale turns out to be much larger than typical sample 
sizes in numerical studies. We find the only coherent object to be a neutral spinon excitation created by the doped hole 
whose energy spectrum is responsible for an isotropic "quasiparticle" dispersion. Thus, our results provide a natural 
explanation of the ARPES results and support the conjecture made in Ref. We also discuss how the observed line 
shapes reflect spin-charge separation, and finally, in the Appendix, how the remnant Fermi surface structure may be 
understood as a peculiar consequence of the phase string effect at higher energies. 



II. PHASE STRING: THE KEY EFFECT INDUCED BY THE MOTION OF A HOLE 



In this section, we examine the effects induced by the motion of a single hole in an antiferromagnetic background. 
In Sec. IIA, we shall review a few basic results of the slave-fermion formalism. We then go on to discuss the effect of 
the phase string on the kinetic energy of the hole, spin-charge separation and the localization of the holon due to the 
phase string. 



A. Slave-fermion formalism of the t J model 



We begin with the slave-fermion representation Qo- — fjbia-{—ay and express the t — J model Ht-j — Ht + Hj 
in the following form: 

Hj = -^Y.(^w^^tj (2.1) 

with 
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(2.2) 



and 



with 



(2.3) 



(2.4) 



At half-filling, the bosonic resonating-valence-bond (RVB) description given by Liang, Doucot, and Anderson |2|] 
has provided by far, the most accurate picture (see Fig. 2) for both long-range as well as short-range antiferromagnetic 
(AF) correlations in two dimensions (2D). In this theory, it is found [^,^ that the long-range AF correlations, 
including the AF long-range order (AFLRO) occurring in the thermodynamic limit, constitute only a small fraction 
of the ground-state energy and the dominant contribution mainly comes from the short-range RVB correlations. 
Thus, the ground state at half-filling may be regarded as AFLRO -|- RVB, of which, AFLRO is the most vulnerable 
part easily removed by either temperature or doping with little energy cost. The mean-field version of this bosonic 
RVB description is known as the Schwinger-boson mean-field theory (SBMFT) |Q. This theory works quite well at 
half-filling, and is characterized by a bosonic RVB order parameter 



A'^ = (Af,) ^ 



(2.5) 



But away from half-filling, the problem is highly nontrivial. The reason can be attri but ed to the fact that the 
hopping integral {Bij) = in TJj. This is generally true due to the sign a appearing in (2^). So the motion of the 
hole will be very sensitive to how the spin backflow Bij is treated which is a non-perturbative problem and is the key 
issue to be dealt with in the present work. 

In SCBA approach the hole can acquire some kinetic energy through the dynamic fluctuations of Bij . In this 

approach, only the long wavelength fluctuations associated with AFLRO is considered, where Bij is approximated in 
large-S* (spin wave) expansion [B| by 



Bi 



B. 



LSW 



= h 



(2.6) 



in which ho denotes the condensed part of the Schwinger boson field, corresponding to AFLRO Now, the hopping 
of the hole is assisted by the fluctuations of B^^w _ xhis is the idea behind the SCBA, and within this approach it has 
been found that the hole behaves just like a Landau quasiparticle, known as the spin-polaron, with four Fermi points 
at ko = (±7r/2,±7r/2) We reemphasize that within this approximation only the long-wavelength fluctuations 

associated with AFLRO are involved, and the dominant short-range RVB correlations, which are independent of 
AFLRO characterized by &o 7^ 0, are completely neglected. 

The validity of the SCBA approach is based on the presumption that spin mismatches, described by the spin 
backflow Bij induced by the hopping of the hole, can be repaired through spin flip process in Hj. Thus, the final 



FIG. 2. The resonating-valence-bond picture [ p2|j2^ of bosonic spins provides a highly accurate description of both 
short-range and long-range antiferromagnetic spin correlations in the Heisenberg model. 
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state of the hole, after traversing a closed path through the antiferromagnetic background, is identical to its initial 
state. This would be true if one focuses only on, say, the z-component of the spins. But since we are dealing with 
a quantum spin system in which each spin has three components that do not commute with one another, there will 
actually be three components in the spin mismatches induced by the motion of the hole. For the spin-1/2 case, it 
has been explicitly proven p^ ] that such string- like spin defects cannot be simultaneously repaired through Hj. Once 
the spin configuration (i.e., the z-component) disordered by the hole hopping is restored by spin flips at low energy, 
a string defect in the transverse components remains, which is described by a sequence of signs in the quantum 
description of the hole. This defect is called the phase string The effect of this phase string, as we shall see in 
the following two subsections, manifests itself both in the expression for the kinetic energy of the hole as well as its 
propagator. 



B. Effect of the phase string on the energy of the hole 

In this section, we shall demonstrate that the total energy at momentum k for the one-hole case can be exactly 
formulated as 

Ek^Eo-^Y.^'^'^'''''''^^^^ , (2.7) 

where 

M,,^ M[c;{^}](-l)^i . (2.8) 

{c}{4'} 

Here, Eq denotes the spin ground-state energy with the hole being fixed at a given lattice site. The energy gain due 
to the hopping arises solely from the second term. In the expression for Mij, the summations run over all the possible 
paths of the hole connecting i and j, {c}, as well as spin configurations, {4>}. The weight functional M is positive 
semi-definite, 

M[c;{(j3}]>0 (2.9) 



such that the phase factor (— l)^-;^ is "uncompensated" in (2^). Here counts the total number of J, spins (or 



I spins by symmetry) being exchanged with the hole moving along the path c. As illustrated in Fig. 3, (— 1)^<^ — 
...(-f 1) X (—1) X (—1) X ... is the phase string on the path c, which was first identified in the exact formulation of the 
single-electron Green's function [Q. Later we shall see that, owing to the factor (— 1)^<: , the lowest energy of E'k is 
obtained for k = kg, consistent with numerical calculations |0-p3. 



From (2/7), (2_^) and (2_^), we come to an important conclusion. Without the phase string factor (— 1)^ , the total 
energy would certainly be lower (as M > 0). Thus, {—1)^^ represents the frustration on the kinetic energy of the 
hole. This effect is (i) irreparable as no other signs can be generated from the spin background to compensate it; (ii) 
singular since a change of iVj by ±1 will lead to a maximal change in (—1)^=^; and (iii) thus expected to dominate 
the low energy dynamics of the hole. A similar effect is well known in fermionic systems where each fermion's path is 
also weighted by a "phase string" determined by the number of fermions "exchanged" with it. The difference in this 
case is that the "exchange" is between two different species, the hole and spins. Therefore, such a phase string effect 
implies an intrinsic mutual statistics between these two degrees of freedom [[L9[ . 



The proof of (2.7)-(p.9|) is straightforward. We use the Wigner-Brillouin formula, 

E^ = Eo + {'i>o{k)\[Ht + HtGjiE^)Ht + ...]'\M'^)) , (2-10) 

where Gj{E) = 1/{E - Hj) and [...]' excludes |-I>o) as the intermediate state. Here |$o(k)) = ^/VNJ2^e''^''''\'^Q^) 
with |$o*'') denoting the ground state of Hj with the hole localized at a site i, viz., Hj\^q ^) = Eq\^q ^). One can 
expand l^o*"*) in terms of the complete spin-hole basis (i))} {4> being a spin configuration with the hole at site i) 

as |$o'^) = '}2<t> X%\4>'i {i))- As shown in Ref. |l^Jl^, x\ ^ 0, which means that the Marshall sign rule still apphes 
to the doped ground state if the hole is not moving. By inserting the complete set of basis states and following the 



steps outlined in mm for calculating the single-hole propagator, one can easily obtain (2.7), with the weight functional 
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K-1 



M[c;{0}] = xUi, n (<^'^'; im,)\PGj{Ei,)P\cj,^; (m,))(-t) , (2.11) 

in which and 0"^^ correspond to spin configurations in intermediate states, and the hole site is on a path c 
connecting two arbitrary sites i and j': mo = i,mi, ...,mK = j {4>^ = 0, 0^ = 4>')- Following Ref. |l^ , we can show 
that {4)"] {ms)\PG j{E)P\(j)^~^ ; {rris)) < as long as E < Eq (note that the projection operator P — 1 - |^'o)(^'o| has 
no effect in the thermodynamic limit N — > oo). Since we are interested in low energy states _Bk < Eq, the weight 



functional M in (2.11) is always positive semi-definite. 



C. Phase string effect and spin-charge separation 

In the above subsection, we showed that the phase string manifests itself as irreparable sign (phase) frustrations 
induced by the hole motion. In this section, we shall establish an intrinsic connection between such a phase string 
and spin-charge separation. 

To this end, let us first consider the origin of the phase string factor (— 1)^<^ in (2.8). Equation (2.4), defining , 
can be rewritten as 

A, = (4)e^^- , (2.12) 

where 
and 

=±(7r/2)[l-a,j] . (2.14) 

Here cr.y = l(^l)i if an | (|) spinon is exchanged with the holon at the link (ij)- For products around consecutive 
links enclosing a loop F, 

^B^,B,u...Bu = (^B%B]^...Bl\^ e'Er^ , (2.15) 

with 




FIG. 3. Phase string as a sequence of signs on the hole path c, determined by the index of each spin exchanged with the hole 
at every step of hopping. 
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j+4>jk+---+4 



(-1) 



where N-^ denotes the total number of | spins exchanged with the hole along the loop T. Thus, the phase 



(2.16) 
of the 



spinon backflow operator (2.12) is the source of the phase string factor in (2 
3. 



This is illustrated by the inset in Fig. 



Since the ground state satisfies the Marshall sign rule plj , it is easy to see that 



n 



> 



(2.17) 



half-filling 



(Note that the Marshall sign is totally gauged away in the definition of bi„ by the sign factor (— cr)* in the slave- 
fermion decomposition (l9j.) This is of course consistent with the previous conclusion that the nontrivial phases in 
the low-energy states all come from the phase string. A hole slowly hopping around the loop F will then acquire a 
Berry's phase <i>r given by 



half-filling 



Imln 



half-filling 



(2.18) 



So the phase string effect generally leads to a nontrivial Berry's phase picked up by the hole, as opposed to the 
quasiparticle picture (of SCBA, for instance) in which the states of the hole before and after traversing a closed path 
are identical. 

Now that we have established the connection between the phase string effect and the spinon backflow, we are in 
a position to discuss spin-charge separation. Let us initially suppose that the hole behaves like a quasiparticle with 
both charge and spin quantum numbers. This means that the holon and spinon should be confined together. In 
this case totally one "f (|) spinon will have to be effectively "transferred" back from the final location of the hole to 
the initial location to ensure a precise spin quantum number | (|) being transported with the hole. It then requires 
that at each step of the holon hopping, the spinon backflow is fully polarized, only to involve a ] (J,) spinon being 
"transferred" backward. Otherwise, since the spinon backflow or the phase string effect cannot be "repaired" , any 
local fluctuations of the spin polarization of the backflow spinous, no matter how weak they are, will be accumulated 
to become arbitrarily large at a sufficiently long path to invalidate that the hole carries a precise spin quantum number 
s = 1/2. In other words, for the spin-charge confinement picture to hold, one must find that = 0, or, = the total 
number of links on any path c, such that the phase string factor (—1)''^=^ becomes trivial no matter how long the 
path is. But based on (2.11) one can easily rule out such a scenario, which imposes the extreme restriction that 
the probability for any other choices of vanishes, since an equal probability for both kinds of spinon backfiow is 
explicitly given in the hopping term ( |2.3| ). In fact, an analysis [^,^ of the weight functional in the single-electron 
Green's function has already lead to the conclusion that the phase string effect must be nontrivial (which actually 
causes the spectral weight Z vanish as shown in Refs. p8| , p^ ). 

Therefore, the irreparable nontrivial phase string effect directly leads to a true spin-charge separation in the single- 
hole doped Mott antiferromagnet. It confirms the conjecture made by Anderson |^]. We note that spin-charge 
separation has also been discussed in some exact diagonalization approaches for momenta close to (0,7r) ||l5| , |l6|l . In 
the following, we will go a step further and discuss some unique features associated with spin-charge separation. 

Let us define the holon propagator Gh in energy space, Gh{E) = (/jG'(i?)//)j^g^j£_gjjjj^g, where G{E) = \/{E — 
Ht-j)- On expanding in terms of Ht, G[E) = Gj + G.jHtGj + we get 



GH{J,^\E)=Y, 




Gj 



half-filling 



(-1) 

half-filling 



(2.19) 



where ^e({fTs}) = (nf=iGj(-t) 



bl 



Gj with A^j^ — Es(-'^~'^s)/2- In (2.1£), c denotes a path connecting 



two sites i and j, and rus in the definition of Sc is a lattice site on any path c. As with the the positive semi-definite 
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FIG. 4. Spin-charge separation in the single-hole case: the injected hole decays into a neutral spinon and a spinless holon, 
and the holon will pick up the full effect of phase string. Representing the spinon backflow effect, the irreparable phase string 
ensures spin-charge separation as discussed in the text. 



functional M in the expression ( p.8|) for £'k, one can easily prove that (Sc) > for energies E < Eq. This means 
that there are no other sources of phases at low energies to repair the phase string factor {—1)^^. It should be noted 
that the single- electron Green's function for the one-hole case has been exactly formulated previously in a very similar 
fashion [Q. There, it was found that each hole path is also modulated by the same phase string (—1)^"=^. This result 
taken in conjunction with ( 2.19| ), shows that the effect of the phase string is entirely in the holon sector, while the 



spinon released from the hole is not directly associated with such phase frustrations. Such a picture is shown in Fig. 
4. This is not entirely surprising, as we do not expect the motion of a single hole to affect the (thermodynamically 
large) spin subsystem, except possibly at short time scales or high energies. 



D. Localization of the holon due to the phase string effect 

Thus far, we have discussed the general properties of the phase string effect in a rigorous form. To further study 
the one-hole problem in a quantitative way, we need a suitable framework for approximation. As mentioned earlier, 
the mean field version of the bosonic RVB theory (SBMFT) provides |Q a fairly good description of the undoped 
antiferromagnet. We now assume, consistent with the discussion so far, that the presence of one hole does not alter the 
low energy properties of the spin subsystem as the irreparable phase string is to be solely picked up by the hole. In the 
following we shall discuss a direct consequence of the phase string effect on the charge part within this approximation: 
the localization of the holon. 



At half-filling, SBMFT characterized by (2.5) is described by the Bogoliubov transformation [24 



bia = -y= ^ ?7ko-e*'^'^''"' Uk7k<T - Wk7k-o- ' (2.20) 
where 7^^ is an operator that creates an elementary spinon excitation with a gapless spectrum at T = 0^, 



= 2.32 1-^2, (2.21) 

and Wk 1/V2{2.32J/E^ + l)i/2, vi, ^ 1/^2(2.32 J/£:^ - l)i/2sgn(Ck), with ^k = (cosfc^^, -I- cosky)/2. Within this 
mean field theory, it is easy to verify that = = 0. If we use this result naively, when there are — 1 

spins (and a hole), we would conclude (erroneously) that i?" = {Bfj) — 0, namely, the amplitude of the hopping 
integral vanishes. This is incorrect because the hopping term connects two spin subspaces, corresponding to the hole 
at site i and j , respectively, which are not necessarily identical in symmetry. This subtlety has to be incorporated in 
any calculation, as will be done below. 

As a first step, let us go back to the case of half filling. Here, it is important to recognize that the factor rj^^ in 



(2.20) which satisfies 77^ — V-a ^-iid \ri\ = 1 cannot be completely determined at the level of mean field theory; i.e., 
the mean field order parameter {bi„hj_„) is independent of the choice of ryko-- Therefore, there is a hidden symmetry 
in SBMFT, which is related to the exact local particle-hole invariance of the system as to be discussed later. 

We now exploit this symmetry in SBMFT for the case when there is one hole. Fixing the position of the hole, one 
may define a subspace for the spins, which, at the level of mean field theory, is again described by the Bogoliubov 



transformation (2.20). It should be remembered that equation (2.20) now describes an A^ — 1 spin subspace. By 



7 



choosing r^kcr = [— S(7n(^k)]'^'' , where fc/j = if the hole is on the even sublattice and kh = 1 ii the hole is on the odd 
sublattice site, one has ryko- — S(7?T.(^k)'7kcr when the hole hops between the sublattices ||2^]. Then noting that in B^^, 
bj^ and bia belong to two different spin subspaces, one obtains 



^° = |El^k|T^E«0.4 . (2.22) 
It is easy to verify that the above choice of rj\^„ optimizes Thus, we obtain a finite amplitude of the hopping 



integral for the motion of the holon between two sublattices. (Note that in the summation of (2.22), k = has been 



removed excluding the bo component. Using (2.4), one can easily see that bo part has no contribution due to the 
sign in a.) The mean field result remains approximately the same for the superexchange term Hj in which the holon 
position is fixed and 77k plays no role. 

This hopping amplitude -Bo 7^ is consistent with A'* 7^ 0, and we may understand this in the following way. 
Physically, each spin subspacc has a hidden local particle- hole symmetry, bia — > ^I-o-- exact and can be preserved 
in the RVB description if the no double occupancy constraint is strictly implemented locally. Let us assume that the 
holon is, say, at an odd sublattice site. Then, on performing a transformation bi^ ^l-a only in the corresponding 
spin subspace, it is easy to see that _Bj?- at the hopping bond transforms into the RVB order parameter A*. Note 
that at the mean-field level in SBMFT, the local hard-core constraint is relaxed where the particle-hole symmetry 
is no longer exact. But it can be checked that the symmetry ryko- — S(7n(^k)'7ko- in SBMFT still approximately 
corresponds to bi„ — > b\_^ at the global level according to (2.2C) (It would be exact if |uk| = |wk|). It tells us that 



the direct hopping term for the holon indeed originates from local RVB spin pairing with a particle-hole symmetry in 
each spin subspace. 



With 7^ 0, Gh in the first line of (2.19) may be written, in the limit oi E ^ —00, as 

G>,{E^-^)^ynl(e'T.J\ (2.23) 



{c} 



half-filling 



where — 2'^(S'c) half-filling' ^o^e that in E ^ —00 limit Gj 1/E becomes commutable with e*"^. We can now 
calculate the phase string average (e'^c"^) half-filling explicitly as 

\ /half-fiUing ^ ' 

^g«ko-(r.-r,)g-i{E^(0-(0»]"> ^ (2.24) 

where " ('^>)]^> = 4 ((E, : - f)^ and : 5f := - (5f). Using the results for {S^S^) from SBMFT M, 



one finds in the large |ry| limit, the asymptotic form for expression (2.24) and thus, the holon propagator as 

G;,(i,i;i?)~e"^«-("----^)e-^ . (2.25) 

The localization length is determined numerically, for the case E = —00, as Xl{E = —00) ~ 2.2a for parallel to 
the x- or y-axis (Here, a is the lattice constant). 



The exponential decay (2.25) of the holon propagator, a typical characterization of localization phenomenon, is 
fundamentally different from the power-law decay for the propagator generally expected when the hole behaves like 
a well-defined quasiparticle. The very fact that the holon propagator behaves like ( ^.25| ) a,s E ^ —00 is enough to 
invalidate conventional quasiparticle behavior. As in conventional localization problems, we expect \l{E) to increase 
with the energy E and to represent the true localization length scale sX E > Eq {Eq is the ground-state energy). 
Later, we shall determine the localization length Xl{E) numerically, in the physical regime of E. We find it to be 
generally larger than the overall scales of sample sizes (< 6a x 6a) used in exact-diagonalization calculations. So, the 
hole localization caused by the phase string effect cannot be directly detected by exact diagonalization, owing to the 
limitation of the sample size. On the contrary, as in conventional localization problems, a localized electron can be 
mistaken for a delocalized electron simply because the sample sizes are smaller than the localization length. Thus, we 
argue that the well defined quasiparticle peak seen in finite size calculations is an artifact of small sample sizes and 
cannot persist at scales beyond Al- 
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III. SPECTRAL FUNCTION 



A. Effective theory 



As pointed out earlier, the bosonic spin RVB pairing gives an extremely good description of the undoped antifer- 
romagnet for both short-range and long-range spin-spin correlations. The direct hopping of the holon originates from 
short-range RVB pairing with a local particle-holc symmetry. Unlike the long-range spin correlations, the short-range 
correlations are not sensitive to doping and the local RVB order parameter A'' provides a certain local "rigidity" that 
underpins the doped antiferromagnet, as long as the the average spacing between the holes is larger than the distance 
between nearest neighbors. Indeed, such a direct hopping term can persist into the metallic (superconducting) phase 
at finite doping, as discussed in Ref. However, unlike the metallic system, the one-hole case is different as we 

do not expect the AF spin background to change at thermodynamic scales. As discussed in the previous section, one 
may, to leading order, assume the spin background to be the same as that at half-filling. The feedback effect on the 
spin part at high-energy, short-distance scales will be considered in the Appendix. 

In the previous section, we saw that the holon picks up the effect of the phase string and that the effective 
Hamiltonian can be written as 

Hh^-thY,e'^''f!fj+h-c. (3.1) 
with th = B^t ~ 0.4t. Here, the nontrivial effect arises solely from the phase (f>ij, reflecting the irreparable phase 



string effect. In terms of ( 2.14 ), we may rewrite 
where 

and ttij = ko • (r^ - r^) [o-y]. 

To characterize the strength of , let us consider the gauge invariant quantity a^, i.e., the fictitious "magnetic" 
flux seen by the holon hopping around a plaquette: i^i + x^i + x + y^i + y^i. One can easily get 

Generally, (X]n '^ij) = 0? ^-iid the strength of the quadratic fluctuations is given by 



(Enajj)^) = TT^ ^1 - ^ (S, • Si+£+j;)j . 



(3.5) 



Using the SBMFT one can estimate ^JdJ^n '^ijY) ~ O.SGvr and find that, for two plaquettes separated far away from 
each other (i.e., the distance R12 ^ a). 



((En«..)i(En«^fc)2)-0(^]^j- (3-6) 
Note that the spatial correlations between the fiuxes threading through different plaquettes fall off rapidly in (^.(j). 



This is because the contribution from the long-range AF fluctuations is strongly canceled out in (3^). Similarly, the 
correlations of fluxes per plaquettes at different times also decay very quickly. Thus, as the first order of approximation 
one may treat as a gauge field describing random fiuxes in the white-noise limit with the strength controlled by ( ^.5| ). 
(The phase string effect beyond the random flux approximation is briefiy discussed in the Appendix in conjunction 
with the momentum dependence of the energy-integrated spectral function). 



In the effective model (3J_), the e ffect from the longer-range spin correlations involving the AFLRO, which infiuence 



the hopping of the holon through (2.6), has been omitted. This process has been considered in the SCBA approach 



as the sole source assisting the holon hopping. In the following we reexamine it in the presence of the bare holon term 
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(|3.l[). With the Bose condensate 60 7^ (AFLRO), one may express cta = (— cr)'cicr + {—crY fj ■ her '■ with Ci„ = 60// 
and : bia- := bi^ — bo. Then, Eq.( |2.6| ) leads to the following hopping Hamiltonian, 

E " + 4 J + '^■c- ■ (3.7) 



If there were no bare hopping (3.1) for the holon, H^^^ would represent a virtual process, which, by the SCBA 



treatment [R-fiB, would give rise to the well-known quasiparticle description for Cia-. But in the presence of ( [3.1| ), 
Hj^^^ literally describes the process for (the "quasiparticle" ) to decay into a holon-spinon pair. In this case, the 
propagator for Cg- will no longer behave like the one for a dressed quasiparticle found in SCBA. In contrast, to 
leading order approximation, it will be simply pro portional to the propagator of : bi„ : to be discussed below. So 
in the following we do not consider H^^^ in (|3.7| ) and simply switch Cia off, by assuming either the sample size to 
be large but finite or T = 0+, such that &o = without loss of generality. 



B. Spectral function 

In this subsection, we shall show that the main features observed in ARPES can be described within our formalism. 
We attribute the broad spectra that are observed to the fact that the physical electron is a convolution of spinon and 
holon excitations. We show that the dispersion of the photoelectron is governed by the dispersion of the spinon, and 
consequently, isotropic. We also show that the phase string plays a crucial role in causing the sharpest spectra at low 
binding energy locating at (±7r/2, ±7r/2) {i.e., the "Fermi points"). 

The single-electron propagator may be expressed in the following decomposition form 

Ge » iGb ■ Gh (3.8) 

where 

Gb{i,r,t) = -i{~ay-= {TtbUt)b]M) . (3-9) 

and 

Gh{l,r.t) = -^{TtI}{t)fM)■ (3-10) 
Here we mainly focus on the properties of the spectral fun ction defined by yl^(k, cj) = — ilmG'e(k, [x))sgn(ti>), which 



can be measured by ARPES. Using the decomposition ( p. 81) , one finds the convolution law 

A^(k, oj) - e{-Lu)^ E f d^'Phi^ - k, C^' - C^)pb(k', LO') (3.11) 
at T = O"*" where pb and ph are the spectral function corresponding to Gb and Gh, respectively. 



By using (2.2C) one can easily obtain 

ImGf, = -Ti[ul5{u - E^) - vl5{u + E'^)]. (3.12) 

As noted above, in the one-hole case this half-filling mean-field propagator should not be affected thermodynamically 
by the motion of the single holon. Then the corresponding electron spectral function is written as 

A^(k, ^) = ^ E ^k'P^Uk' - k, -E^, - u) (3.13) 

k' 

where the condition ph 7^ only at w > is used which defines the bottom of the holon band at a; = 0. 

The overall feature of ^"^(k, w) versus uo is shown in Fig. 5 at different k's along (0, 0) to ko. The whole energy 
range of the spectral function is about 16J 2.2AeV (for J = Q.lAeV). Here we have chosen = 2J {t — 5J) in 



(3.1) in obtaining the holon spectral function p^ which is calculated numerically by exactly diagonalizing under 
random fiux with | < 0.867r in the white noise limit. The sample size is 32 x 32. 

Even though not delta-function-like, the spectral function do show peaks which become sharper near both energy 



bottom and top as the momentum approaches kg. According to (3.1), the bottom and top of the holon band are near 
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FIG. 5. The full shape of the spectral function calculated at k = (n /2,ii /2), (7r/4, 7r/4), and (0,0). Note the whole energy 
range extends over 16 J « 2.2AeV (for t = 5J, J = Q.UeV). 



the momenta shifted away from (0,0) and (tTjTt), respectively, by kg, where the spectral function ph is sharpest in 
contrast to the broadest feature near the band center due to the random flux. On the other hand, the spinon spectrum 
is maximum at kg and vanishes at momenta (0,0) and (tt, tt). The convolution law of (3.11) then combines the 
contributions from the holon and spinon to result in Fig. 5. Note that the peaks are asymmetric at the top and 
bottom of the band. The effect of the spinon spectrum is most visible at lower binding energies, as will be further 
discussed below. A similar asymmetric structure in one-dimension has been found in Ref. [ p8[ . 

Let us now focus on the structure of A'^ at low binding energies in comparison with ARPES measurements on 
Sr2Cu02Cl2 ||l|,|l and Ca2Cu02Cl2 In the left panel of Fig. 6, we plot A"" within the energy range of 8 J ~ lAeV 
at k positions along a Une from (0,0) to {tt/2, tt/2) and then from {tt/2, tt/2) to (tt, 0). The peak or edge position 
shows an isotropic dispersion along (0,0) to {n/2, tt/2) and from (7r/2, 7r/2) to (tt, 0), consistent with the experiments. 
Such a dispersion has a bandwidth about 2 J and is clearly correlated with the spinon spectrum ii'k+ko- "^^^ latter is 
marked by small bars at different k's where the position of the minimum -B^+ko fixed around lo = 0.7 J as a 

reference point. In the right panel of Fig. 6, a plot of A'^ along (0,0) and (0,7r) is shown where the peak is replaced 
by an edge which looks dispersionless, also in good agreement with the experiment. Indeed, a very flat (only about 
13% change) E^+i^g along (0,0) to (0, tt) is indicated in the figure by the small bars. In Fig. 7, the spinon spectrum 
^k-i-ko fits observed ARPES "quasiparticle" spectrum data reasonably well over the whole Brillouin zone along 
ko (S) to (0,0) (F) and (0,7r) {X) (which are symmetric in our theory) with J = 0.14ey. Note that the low-energy 
scale is determined by J instead of t and the features shown in Fig. 6 are not sensitive to t. Thus, second - or further 
neighbor hopping terms are not expected to play as important a role as they do in determining the energy bottom, 
Ek, shown in Fig. 1. 

The reason that the low energy peak or edge of the spectral function is correlated with the spinon spectrum can be 
easily understood by notin g the fact that ph(k,uj) becomes the sharpest near the bottom = with k sa ko which 
contributes to A^ in ( 3.13| ) at k' « k + ko and cu ~ — E'k-fko- In Fig- 8(a)j the spectral function along (0,0) - ko 
is shown when the maximal strength of | | is reduced from O.SGtt to O.Itt. Here the correlation between the 

low-energy peak or edge positions and the spinon spectrum is even more apparent as the holon spectral function ph 
becomes sharper near lu — and ko. The evolution from a peak to an edge as k moves away from ko is because the 
spectral function of the holon is quickly broadened away from the band bottom in the presence of the random flux. 
Note that compared to the experiment the peak structure in Fig. 8(a) is a bit too sharp which implies that the actual 
strength of | aij\ should be closer to our estimation used in Fig. 6. 

The line shape of the spectral function shown in Fig. 6 and Fig. 8(a) looks strikingly similar to the ARPES 
measurements in Sr2Cu02Cl2 and Ca2Cu02Cl2 |^. It reflects essentially the convolution law of spinon and holon, 
i.e., the spin-charge separation. But we would like to point out that the simple convolution law is not enough. Here 
the coherence of the spinon in contrast to the incoherent holon is key to the characteristic features of the line shape 
and the "dispersion" of the peak or edge structure. To illustrate this point, we replace the coherent factor in ( 3.13| ) 
by a constant (i.e., v^. — 1), and re-calculate the spectral function which is plotted in Fig. 8(b) along (0,0) - ko. One 
sees that the whole low-energy peak-ed ge fe ature in the left panel of Fig. 6 and Fig. 8(a) totally disappears in Fig. 
8(b). In this case, the summation in ( 3.13 ) smears out the peak structure and A'^ is more or less like the density 
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FIG. 6. The spectral function at low binding energy along different momentum scans. The small bars mark the spinon 
spectrum -Bk+ko ~ Mo with fio = 0.7 J. 



of states for the holon. So the coherent factor in the spinon propagator is crucial for the "quasiparticle" peak to 
emerge in the spin-charge separation formulation of the spectral function, which is the unique result of the bosonic 
RVB description of spins. This explains why the correct line shape has not been directly obtained in the slave-boson 
formalism, even though the prediction that the "quasiparticle" spectrum obtained in ARPES actually corresponds to 
the spinon spectrum was first made 1^ there. 

The incoherence of the holon is also an important factor as discussed above. It is well-known that a particl e mu st be 
generally localized, at least in the low-energy regime, in the presence of the random flux governed by Hh in ( |3.l| ). We 
have checked the localization lengths corresponding to the case of Fig. 6 using standard methods ||2^ . The localization 
length quickly jumps from 4-15 (in lattice units) to 36 — 50 near the holon band edge which covers an energy range 
related to the low-energy peaks of the spectral function in Fig. 6. As noted before, the localization length scale 
of the holon is generally much larger than the maximum sample sizes (< 6 x 6) in numerical exact-diagonalization 
calculations |^2|. Thus, numerical calculations providing a "metallic" quasiparticle picture should be irrelevant to the 
true picture at large length scales. 

Finally, to conclude our discussion of ARPES in the Mott insulator, we turn to the experimental results of Ronning 
at al. These authors identify a "remnant Fermi surface" in the insulator from ARPES measurements: they define 

700 
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FIG. 7. The "quasiparticle" spectrum determined by the ARPES |^ (open square) and by the spinon spectrum 
(solid curve). 
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FIG. 8. (a) The spectral function similar to the left panel of Fig. 6, but with the random flux strength per plaquette reduced 
to O.Itt (see text). Small bars marked the same spinon spectrum as in Fig. 6 (but with fio = 0.2J); (b) The line shape of the 
spectral function in the left panel of Fig. 6 is greatly changed if the spinon coherent factor is replaced by a constant 1. 

the so-called relative momentum distribution 

"fc(^o)= T— A^(k,a;) , (3.14) 

with the cutoff energy ojq = —0.5eV, then identify "kp" by locating the position where drops suddenly. Then a 
"remnant Fermi surface" is found as the contour of steepest descent of n]^. Such a "remnant Fermi surface" roughly 
coincides with the large Fermi surface expected for a free electron gas at half-filling concentration. It should be 
emphasized that this contour is not an equal energy contour. In fact, the peak (or edge) of the spectral function 
disperses at momenta along the "remnant Fermi surface" as much as that shown in Fig. 7. 

So the "remnant Fermi surface" mainly indicates a strong momentum dependence of the energy-integrated spectral 
function, n^. We believe this is due to the effect of the mobile hole on the spin background (at very short time 
scales). So far we have not considered the influence of the phase string effect on the spin degrees of freedom. This 
is based on the fact that a single hole does not affect the spin background at the thermodynamic scale. However, 
the AF background surrounding the doped hole should be still strongly distorted by the hopping at i ^ J, which 
in turn will feed back on the hole itself through the phase string effect. Obviously, this happens at short length and 
time-scales. Consequently, one has to go beyond the effective Hamiltonian ( |3.1[ ) as well as the random flux treatment 
of the phase string effect. Not being within the main purview of this paper, we relegate the discussion of such effects 
to the Appendix. 

IV. CONCLUSIONS AND DISCUSSIONS 

Although antiferromagnetism in the Mott insulator has been well understood based on the bosonic RVB description, 
the problem of a single hole doped into such a system is found to be nontrivial. Unlike the spin-polaron picture in 
which the hole is predicted to behave like a Landau- Fermi quasiparticle carrying a flnite-size spin distortion around it, 
the real picture for the motion of the doped hole is that it always picks up a string of signs from the spin background 
during its hopping, known as the phase string, which leads to the destruction of the "coherence" of the doped hole as 
a Landau-Fermi quasiparticle. 

The phase string characterizes the effect of spin backflow during the motion of the hole. The irreparable nature of 
phase string is intimately related to spin-charge separation: the hole cannot carry a precise spin-1/2 quantum number 
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with it. In contrast, if the holon and spinon were tightly confined together, the phase string effect would become trivial 
(reparable), as is the case when one naively uses the c-operator to replace the spinon-holon decomposition in ( ^.5| ). 
Hence, the demonstration that the phase string effect cannot be "healed" through spin superexchange interaction 
at low energies, presented in Sec. IIB and Refs. points directly to spin-charge separation. It is also found 

that the spin backflow, or equivalently, the phase string leads to mutual statistics between the hole and spin degrees 
of freedom. That the backflow of spin current could lead to mutual statistics and spin-charge separation was also 
anticipated by Baskaran . 

In the one- hole case, the phase string effect mainly acts on the holon part, causing its localization. But the spinon 
remains coherent in the bosonic RVB background. Consequently, during its propagation, a bare hole releases the 
coherent spinon constituent whose isotropic dispersion essentially determines the position of the low-energy peak or 
edge of the spectral function. On the other hand, spin-charge separation characterized by the convolution law of the 
holon and spinon propagators is responsible for the intrinsic broad feature of the spectral function at higher energy. 
In the phase string formalism, the singular part of the phase string effect also decides the Fermi points ko in the 
low-energy, long-time limit as well as the large "remnant Fermi surface" structure in the high-energy, equal-time 
limit. We thus obtain a systematic and natural explanation for the ARPES measurements within the framework of 
the t — J model. Conversely, one may say that the ARPES experiments have clearly presented evidence for spin-charge 
separation, as emphasized by Anderson. 

We find that the "coherence" of the spinon, reflected by the factor in SBMFT, plays an important role in 
determining the ARPES line shapes. A naive convolution without this factor is not enough to explain the line shape 
observed by ARPES. We also note that the strong temperature dependence of the ARPES line shapes over a wide 
energy range [0, may be easily understood based on the temperature dependence of v^. Details will be presented 
elsewhere. Finally, we point out that the temperature dependence of is also crucial to get the correct behavior 
of the relaxation rate (1/Ti) in nuclear magnetic resonance measurements, spin-spin correlation length, and other 
magnetic properties at half-filling within the SBMFT. 

The prediction that the holon as the charge carrier is localized by the phase string effect is also consistent with 
the cuprate superconductors. It is known that the phase with AFLRO (weakly doped regime) is always an insulator. 
According to our theory, the localization is due to the intrinsic nature of the doped Mott insulator instead of external 
reasons like the Anderson localization in the presence of impurities. 

Note that the localization of the hole is because the phase string is solely picked up by the hole, and the antifer- 
romagnetic background remains the same as at half-filling, thermodynamically. At finite density of holes, both spin 
and charge degrees of freedom will be affected by the phase string, which can lead to a metallic (superconducting) 
phase 1 27 1 without AFLRO. We emphasize that the irreparable phase string effect always exists, even in the metallic 
phase. This is because the phase string reflects the competition between the hopping and superexchange interactions 
in the t — J model, which has nothing to do with the existence of AFLRO. Thus , the p hase string effect will persist 
at any flnite doping so long as short range antiferromagnetic correlations persist |l9| , |27|l . 
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APPENDIX A: SHORT TIME SCALE EFFECTS AND THE "REMNANT FERMI SURFACE" 

The holon is localized by phase string in space. But at a short scale within the localization length, the holon may 
behave like in the metallic phase at finite doping where holons are mobile. It has been found |l^,|2^ that in the 
metallic phase the phase string effect will influence both the spin and charge degrees of freedom in such a way that 
the correct decomposition form for the electron becomes 

c,^ ^ hlk^e'^''""\ (Al) 

where h] and bier are the "true" holon and spinon operators and the phase shift field e'^*^ " precisely keeps track of 
the non-repairable phase string effect. In the following we will use this formalism to describe the local, high-energy 
regime in the one-hole case and show a nontrivial consequence of the phase string effect. 
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In the one-hole case, the phase shift field can be written as e 



= (— a-)*e2*i in which 



<i>^-^Ini In {z,-zi) [Y,an\^~l\ 



(A2) 



with Zi = Xi + iyi representing the complex coordinate of a lattice site i and n^^ denoting the spinon number at site I. 

In order to understand n^, let us first take o^o to — oo. In this high-energy or equal-time limit, n^. reduces to the 
momentum distribution nk which is given by 



(A3) 



in which the following expression for the equal-time phase-string factor is used: 



5i5*.(*)e-*5*'(0) 



with 



p Jr 



Z X (r — r;) 
|r - r;|2 



(A4) 



(A5) 



where the path F connects two lattice sites i and j without crossing other lattice sites. Clearly, in the equal-time 
limit, the momentum structure will be mainly decided by the oscillation of the phase string factor given in ( A4). 

J 



For the one-hole case, nk in (A2) is actually trivial and featureless {— 1/2). This is because (cJ^Cjo-) — 5ij{b\^b, 



using the non-doublc-occupancy constraint in which no propagation of the holon at a finite distance can occur at 
strictly equal time on the half-filling ground state. On the other hand, nJJ for a finite cutoff ojo will involve a finite-time 
holon propagation over some finite distance. In this case, the phase string factor in (A4) will show up to determine 
the basic feature while the rest of the propagator involving holon and spinon will mainly give rise to a broadening 
in the momentum space depending on how far the holon and spinon constituents can travel under the cutoff loq. At 
LjQ — — oo, the momentum broadening would reach infinity (as only i = j contributes) such that any momentum 

structure arising from e'-^r'''" '^ gets smeared out. But at a large but finite \ujo\ the momentum structure due to 
the phase string factor should show up. Here as an approximation the equal-time phase string factor is still used so 
long as I Wo I is sufficiently large. 

In the spin channel, a line-integral expression similar to (A4) also appears in the spin-spin correlation function, 
leading to the so-called incommensurate antiferromagnetic peaks. Here we can borrow the similar method used in 
Ref. [^2| to manipulate the contribution of the phase string factor (A4) in n^, which gives rise to four "incommensurate 
peaks" at (zLttk, 0) and (0, ±ttk) in momentum space with k ~ 1 |33|. If this oscillating factor solely decides the 
momentum structure of nj. with the rest term in the propagator mainly contributing a broadening as discussed above, 
then one gets a contour plot of in Fig. 9 in terms of the superposition of aforementioned four peaks (with k — 0.75) 
at an arbitrary broadening (which is presumably controlled by the energy cutoff loq: when loq —oo, the broadening 
in the momentum space should go to infinity such that rik — 1/2). Note that the values of here are only meant 
for relative comparison. 

We see that the overall feature is in qualitative agreement with the experimental data 



whose one quarter portion 



is also shown in the inset of Fig. 9 for comparison. Here we remark that experimental data 



of nl seem 



to be photon-energy dependent, indicating the effect of the electron-photon matrix element, and there also exits 
controversies [^,^ on whether the "remnant Fermi surface" defined near the sharpest drop of in k-space is 
meaningful or not. But in the present approach the momentum dependence of does not represent any real Fermi 
surface structure of the Mott antiferromagnet at half-filling. It only refiects the superposition of four peaks near 
(±77, 0) and (0, ±7r), an enhancement entirely coming from the dynamic effect of the doped hole as the result of a 
careful handling of the phase string effect in the local, high-energy regime. As a prediction, if experimentally \ujo\ is 
taken to sufficiently higher energy, the momentum dependence of should get weaker and weaker and eventually 
—f Uk ^ 1/2 at ujQ — )■ —oo. 
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FIG. 9. The contour plot of the electron momentum distribution exhibiting a "remnant Fermi surface" structure. The 
experimental data M are shown in the insert at the upper right corner. 
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